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1.0 Introduction 


It is estimated that more than 500 x 10® metric tons of dust are 
transported annually by the wind (Peterson and Junge, 1971). Dust storms 
reduce visibility on highways and are responsible for loss of life and 
property through many accidents each year. Atmospheric dust, whether raised 
by winds or injected into the atmosphere by volcanic processes, can also have 
a significant effect on temperature. Thus, windborne particles can have a 
direct effect on the climate. In addition, windblo^gn sands cause abrasion and 
erosion of natural and manufactured objects and encroach upon cultivated 
areas, turning productive land to deeert — a process termed desertification. 

The problem of desertification is enormous and is recognized on all inhabited 
continents of Earth. Agricultural land damaged by wind erosion in the United 
States alone varies from 400 to 6,000 km^ per year (Kimberlin et al., 1977). 

The key to understanding problems associated with windblown sediments, 
including those aspects of desertification dealing with wind processes, is 
knowledge of the physics of windblown particles. 

Many paramet.ers, such as grain size and wind speed, must be considered in 
assessing the entrainment and transport of windblown grains. Beginning with 
Bagnold's (1941) research in the 1930's, numerous Investigators have used wind 
tunnels to analyze particle entrainment. Wind tunnels have the advantage that 
individual parameters can be closely controlled and isolated for study and 
analysis. Gravity is one of the most critical parameters in the analysis of 
windblown sediment movement; however, there is no effective means for 
isolating and assessing its effect during experiments conducted on Earth. 
Hence, a space station would afford an opportunity to conduct experiments in 
zero-gravity and variable-gravity environments. Such experiments would enable 
not only this parameter to be assessed, but through its elimination as a 
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factor in the experimenUB , other critical parameters could be evaluated as 
well. 

The objective of this feasibility study was to determine if a suitable 
apparatus could be designed to analyze aaolian processes for operation in 
space and to assess the feasibility of conducting meaningful experiments to 
address, key (aspects of aeolian processes. To meet this objective a prototype 
apparatus was fabricated and some limited experiments were run to determine 
its suitability for this application, 

^ ^ Pobsntial Experiments in the Spacfe Station Environment 

We envisage at least three general types of experiments that could be 
carried out aboard a space station; threshold studies, swirl (dust devil) 
experiments, and analyses of windblown particle trajectories. In this section 
we consider how experiments in a zero~g environment could advance knowledge of 
aeolian processes. 

2 , 1 Threshold Experiments 

Threshold (u ) defines the minimum winds required to initiate particle 
*t 

motion and is the fundamental factor in aeolian processes. In the 
determination of a general expression of the threshold wind speeds for small 
(~ sub-millimeter) particles, the effect of aerodynamic forces tending to 
dislodge a particle from a bed of loose granular material is equated to the 
effect of opposing forces, namely the particle w<=lght (W) and interparticle 
cohesion (I). The relative magnitudes of these forces have been deduced only 
approximately from wind tunnel tests of threshold speed (Iversen et al. , 1976; 
Greeley et ai. , 1980a; Iversen and White, 1982). 

The elimination of particle weight in the threshold force equation — as 
could be accomplished by conducting experiments in a weightless environment — 
would enable a more accurate assessment of the other factors. The equation of 


Table 1. Nomenclature 


pressure differences 
proportionality coefficient 

drag coefficient 

lift coefficient 

surface stress coefficient 
for a vortex 

aerodynamic drag 

particle diameter 

gravity 

cohesive force 

central laye?: speed 
coefficient (Eq. 8) 

lift 

moment 

particle friction Reynolds 
number 

inner drum radius 

outer drum radius 

radial coordinate, carrousel 
wind tunnel 

characteristic vortex radius 

radius at boundary between 
central layer and outer 
boundary layer 


radius at boundary between 
central layer and inner 
•boundary layer 

U wind speed 

u inner drum surface 

friction speed 

u outer drum surface 

friction speed 

u threshold friction speed 

*t 

Vj. wind speed relative to 
particle 

W particle weight 

0 angular coordinate 

V kinematic viscosity 

p atmospheric density 

pp particle density 

oj tensile stress due to 
interparticle cohesion 

u inner drum rotational 

speed 

Wq rotational angular 

velocity of a vortex 
generator 


equilibrium for a email particle at threshold is 

Da + Lb + M - Ic + Wb 


Pi':.'; :j 

OF 




( 1 ) 


VHiere D, L, and M are aerodynamic drag, lift, and moment, respectively, W is 
particle weight, I is cohesive force (see Table 1), and a, b, and c are 
distances from lines of action of the forces to the overturning point. Figure 1 
illustrates the force system acting on the particle at rest.. Elimination of the 
weight term would aid in the determination of the form and magnitude of the 
cohesive force term at the moment of threshold. 

The equilibrium of the forces at incipient motion can be used to express an 
equation for threshold. For values of particle friction Reynolds number, R^, 
between 0.3 and 10 (normal range for sand on Earth), the resulting 
semi-empirical equation (Iversen and White, 1982) is 
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P„gD 

0.0166 -H— E. 
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1 + 


0.006 

2.5 
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0.092 

(1.928R* - 1] 


( 2 ) 


With the elimination of the weight term, the equation becomes 
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u. 


*t 


0.092 1.5 

0.0001 [1.928R. - 1] T pD 

* p 


(3) 


Experimental data for threshold at one Earth atmosphere are shown in Figures 2 
and 3. The expectation of threshold at zero-g is shown in Figure 4. 

2.2 Swirl Threshold Experiments 

Cyclonic atmospheric motions, or "dust devils", are effective in raising 
large quantities of very fine-grained particles in local settings such as 
cultivated fields. Such motion often occurs in otherwise calm atmospheric 
conditions and may be responsible for the onset of some major dust storms. The 
physics of particle entrainment by dust devils is virtually unknown, although 
some preliminary small-scale experiments have been conducted previously (Greeley 
et al. , 1981). A seml-empirlcal equation was derived from these experiments 
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where is a representation of maximum tangential speed within the 

t 

dust-devil vortex. If the weight force is eliminated as could be accomplished 
in a zero-g environment, the result is 


(r u ) 
o 0 


pCflp + 


(5) 


Elimination of the weight as a force simplifies the equation and would 
allow more accurate measurements of aerodynamic and interparticle forces. The 
effect of the difference between Equations 4 and 5 is illustrated in Figure 5. 
2 . 3 Particle Trajectories 

Once particles are airborne, the forces ,?,cting on individual grains in 
their trajectories are the particle weight and the aerodynamic lift and drag. 
The two-dimensional equations of motion of a particle moving through the air 
can be written as 


4 (3 

D X 

3 p p 


yV C - (x - u) V C 
r L r D 


4 P„ . 

■s — E. D y = -(x - u)V C. - yV C_ 
3p p-^ “^rL •'rD 

2 . 2.2 

*= (x - u) + y 
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( 6 ) 


The dotted line surrounds the weight factor. Experimental and calculated 
trajectories for one-g conditions are illustrated in Figure 6. With the 
elimination of the weight factor under zero-g, the only forces remaining are 
aerodynamic. Thus, experiments conducted in zero-g would enable aerodynamic 
lift and drag to be assessed. 
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^ ® Carrousel Wind Tunnel (CWT) Design 

Boundary-layer wind tunnels appropriate for Investigations of aeollan 

processes typically exceed 10 m In length (e»g. Bagnold, 1941; Greeley et al., 

1981). Because of the probable size constraints In an orbital laboratory, we 

» 

evaluated various wind tunnel designs that would provide a compact apparatus 
and which would be suitable for coLiductlng the types of experiments outlined 
In Section 2, Based on these considerations, we designed a wind tunnel not 
previously u?;ed In aeollan studies. The type of wind tunnel we propose 
consists of two concentric rotating drums (Fig. 7). The space between the two 
drums comprises the wind tunnel test section. Differential rates of rotation 
of the two drums would provide a wind velocity with respect to either drum 
surface. Rotation of the outer drum provides a "pseudo” gravity ("pseudo" in 
the sense that a gravity force acts on the particle only when it is resting on 
the outer drum surface). 

In order to test the concept of this wind tunnel design, a 1/3 F,cale 
model Carrousel Wind Tunnel (CWT) was constructed and calibrated (Figs. 8-10). 
In this prototype, only the Inner drum rotates, whereas in the final 
configuration, bof■^ drums would rotate. The drum walls are plexiglass to 
allow observation of the test section. The inner drum is 712 mm in diameter 
and 258 mm wide and is spun around a horizontal axis at a constant speed of 
i‘15 rpm (53.9 rad/ 2 ), which gives a drum-surface speed of 19.2 m/s (63 ft/s). 
The outer drum is 1062 iiam in diameter and 263 mm wide. In a full-scale model, 
both drums would rotate at controllable, variable speeds. The outer drum is 
se,aled along its periphery, but there is a small gap between the sides of the 
inner drum and the outer drum. It is assumed that a zero-g environment would 
be available on board a space station, and that atmospheric pressure 
environments over the range of 1 bar down to low (~1 mb) pressure would be 
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available. 
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4.0 Prototype Testing and Experimental Comparisons 

A series of experiments was conducted to determine the boundary layer 
properties of the prototype CWT. Taylor (1935) hypothesized nearly potential 
(inviscid) flow for the turbulent flow between concentric rotating cylinders 
and noted that boundary layers near the outer and inner walls should be 
govsrrjed by Prandtl’s mixing-length theory. In CWT it is important that the 
mixing-length theory govern the boundary- layer flow adjacent to the curved 
cylinder wall surfaces because the same theory governs the flow adjacent to a 
plane surface and would be comparable to natural conditions and to conditions 
used in previous threshold experiments (Greeley et al. , 1976, 1980b; Iversen et 
al. , 1976). Experiments were performed in CWT to ascertain if these assumptions 
are correct and if Taylor’s hypothesis is valid. For cases in which only the 
inner cylinder rotates (as in the prototype CWT) and assuming that the surfaces 
of the cylinder walls are aerodynami cally smooth, the following equations for 
the flow between two infinitely long cylinders (Fig. 11) can be derived: 

inner layer (Prandtl boundary layer) 

U »= -• u^^ {2,5 In l(r - F.j|^)u^^/v] + 5»5} 

for R. /R + (v/u*.R )e°*^’^i“^“*l (7) 

10 ”10 


central layer (potential inviscid layer) 

U « KR.w(R + R,)/4r 
1 o i 

for r^/R < r/R < r,/R 
2' o 0 1 0 

outer layer (Prandtl boundary layer) 

U . {5.5 + 2.5 In 1(1 - |- > \ V"” 

0 

for r,/R < r/R 

1 p 0 

< 1 - 0.1108/(R^ 


( 8 ) 


( 9 ) 
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To determine the ^low characteristics in the CWT, wind velocity profiles 
were obtained in the prototype using a TSI model 1010 hot-wire probe anemometer 
at four positions in the tunnel; 17.5, 39.5, 77,5, and 126,5 mm from the side 
of the drum and within about 5 mm from the inner drum to about 5 mra from the 
outer drum. The velocities and distances were non-dimens ionalized by dividing 
the wind velocity by the inner drum rim speed (U/Rj^w) and the distance by the 
outer drum radius (t/Rq). 

The theoretical Eqs, 7,8,9 were then matched to the experimental data by 

first fitting the data in the central layer to Eq, 8. Then the values of r^ , 

r„ , u , u were calculated by trial and error by matching values of speed U 
^ *i *0 

and slope dU/dr at the inner layer - central layer and outer layer - central 
layer boundaries. The results are shown in Table 2 and Figures 11-14. 



Table 2. Speed Profile 

Characteristics 



Station No. 
(mm) 

r^/Ro 

Tg/Ro 

u 

*i 

( cm/s) 

u 

*0 

( cm/s) 

K 

U 

*RiW 

Ri + ^ 
2 

17.5 

0.9016 

0.7915 

46.9 

29.4 

0.316 


0.378 

39.5 

0.9015 

0.7941 

47.2 

29.2 

0.313 


0.37'i 

77.5 

0.9015 

0.7942 

47.2 

29.1 

0.313 


0.375 

126.5 

0.9016 

0.7911 

46.9 

29.4 

0.317 


0.379 


Both the data in the table and in the figures show that the flow is quite 
uniform across the breadth of the wind tunnel cross-section. Figure 15 is a 
comparison of theoretical velocity distributions. In Figures 11 and 13, the 
central and outer theoretical layers give remarkably close fits to the data, 
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and in Figure 14 the central and inner theoretical profile layers fit the data 
quite well. The primary reason for the slight dlscrepancj' between theory and 
experiment, especially for the inner layer in Figures 11-13, is the existence 
of secondary flows in the wind tunnel cross-section which are not taken into 
account in the theory. 

These experiments show that Taylor’s hypothesis appears to be correct. The 
results also show that surface friction spoeds can be predicted a priori with 
some degree of accuracy, facilitating the design rf CWT and enabling a 
reasonable experimental matrix to be developed. 

The velocity profile data of Figures 11-15 show that the desired 
two-dimensionality of the flow in the carrousel is very nearly realized over 
most of the cro8S-sectlon<, The variation of measured speed across the 
cross-section width is Illustrated in Figure 16 for three radial positions and 
shows that secondary flows are concentrated near the inner (rotating) cylinder . 
Turbulence levels of 6Z to 10% were raeaeured within the inner and outer 
boundaries, which are within levels to be expected in normal, flat-plate 
boundary layers developed within conventional atmospheric wind tunnels. 

From these calibration experiments and flow-field determinations, we 
conclude that the CWT boundary layer environment is appropriate for experiments 
dealing with threshold and particle trajectories. 

^ ® Parti cle-Traj e ctory Calculations 

In order to assess the suitability of CWT in the analysis of the 
trajectories of windblown particles, a series of calculations was conducted 
comparing cases for one-g with those of zero-g. The equations of motion for an 
airborne particle, assuming no lift force, are (in a polar coordinate system): 
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V 




» • 


•2 3pC«r 

r-re - gcosQ-h (-rr-f-) V •• 0 

P P 

•< • • 3pC_ • 

r 6 + 2 r 0 + g sin 6 - ■ | > - ) [U(r) - r 0] V^. « 0 (10) 

'^P P 

• 2 * 21/2 
Vp » {r + [U(r) - r 0J } 

Equations 10 were solved for several example cases. • The drag coefficient, Cjj, 
is a function of Reynolds number, assuming a spherical particle (VHilte et al. , 
1975). Figure 17 illustrates particle trajectories in CWT for a zero-g 
vacuum for 5 different initial radial velocities. The coordinate system is 
fixed to the particle launch point and rotates with the outer cylinder. In 
inertial space the trajectories are straight lines, but relative to an 
observer standing on the launch point of the rotating outer drum, as plotted, 
the trajectories are curved. For comparison, a particle trajectory is shown 
in Figure 18, also for zero-g but in an environment of one atmosphere air. 

The initial inward radial velocity of the particle is assumed to be equal to 
the surface friction speed of the outer cylinder (32.6 cm/s). The assumed 
wind speed profile for the calculation is shown in Figure 19. The 
trajectories in Figures 20 and 21 were calculated in order to compare 
trajectories with and without the effect of gravity. The same velocity 
profile as in Figure 19 was used to determine the particle trajectory on the 
Earth’s surface and is shown in Figure 20. Also shown in Figure 20 is the 
trajectory of Figure 18, plotted in a rectangular coordinate system. The 
effect of gravity is to lower the trajectory height by a factor of about 
one-thiri. Since the only force acting on the particle in CWT is aerodynamic, 
the significant difference betv;een the two trajectories should enable much 
more accurate determination of the aerodynamic forces (drag and lift) than is 
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possible in an Earth-based facility. This difference can also be shown ly 
calculating the particle trajectory in the carrousel when it is positioned on 
the Earth’s (surface, as in the prototype tests. For the trajectory of Figure 
21, the particle is initially located at the lowest position of the outer 
cylinder and gravitational attraction is in the negative vertical direction. 
Again the trajectory peak (3 mm) is about one-third lower than that in Figure 
18 (5 mm). 

We conclude that CWT could yield significant data on the trajectories of 
windblown particles. 

6.0 Summary and Conclusions 

A wind tunnel apparatus was designed as appropriate for use in an orbital 
space station. A 1/3 scale prototype system was fabricated and its flow 
characteristics were assessed. Preliminary results show uniform flow and 
boundary layer properties that are in agreement with theory. Experiments were 
conducted in the prototype to determine the feasibility of studying various 
aeolian processes and the results were compared with various numerical 
analyses. Several types of experiments appear to be feasible utilizing the 
proposed apparatus: 

Threshold Experiments under Varying Gravity 

Aeolian threshold experiments for a wide range of atmospheric pressures 
(0.004 to 30 atmospheres) have been conducted previously (Iversen et al,, 1976,* 
Greeley et al., 1976, 1980a, b). All of these experiments have been conducted 
under conditions of Earth’s gravity. It would be extremely valuable to extend 
the matrix of experiments to include values of artificial gravity above and 
below that experienced on Earth, This could be accomplished in CWT by placing 
a bed of particles on the inner surface of the outer drum and rotating the 
outer drum at different rotational speeds. The rotating outer drum provides an 


acceleration directed radially outwards, normal to the surface, thus creating 
artificial gravity. While rotating the outer drum at a constant rate to 
maintain a constant value of artificial gravity, the inner drum speed can be 
changed to increase the value of outer-drum wind-friction speed until the top 
layer of particles leaves the surface at which point the threshold wind 
friction speed can be ascertained. 

Threshold Experiments for Zero Gravity Conditions 

Zero-gravity threshold experiments are valuable because of the elimination 
of the weight term in Equation (1). Such experiments would enable cohesion 
forces to be assessed, as discussed in section 2. These experiments would be 
conducted by rotating only the inner drum, accelerating its speed until 
uticeshold is reached. 

Swirl Threshold Experiments 

Threshold experiments and particle motion in the zero-g environment would 
also he out in a dust-devil simulator, an apparatus used in previous 

experimenus (Greeley et al. , 1981). The same advantages are found in these 
experiments as for threshold experiments, i.e., forces of cohesion and 
atmospheric lift and drag can be separated from the effects of particle weight 
and thus can be determined much more accurately. 

Experiments to Study Particle Trajectories 

Analysis of particle trajectories in a zero-g environment would enable the 
determination of the aerodynamic forces on windblown particles using high-speed 
motion picture obtained during the experiments. The lift and drag forces would 
be determined by measuring particle accelerations, particle speeds, and wind 
speeds, and applying Equation 6 to the results. 

Analyses of Bed-Forms 

Small bedforms, such as ripples, are of great interest in aeolian 
phenomena. Because the physics of ripple formation are poorly understood, it 


is not known if ripples would form in the space station carrousel, although it 
is expected that ripples would form because the particle trajectories relative 
to the bed would be similar to those on Earth (Fig. 20), 

By rotating the inner drum at a faster rotational speed than the outer 
drum, a wind veli;»dLty profile could be generated which would cause saltation 
from a bed of particles situated on the inner surface of the outer drum. By 
observation under a variety of rotational rates and rate differences, 
experiments in the controlled environment of the carrousel could increase 
understanding of the physics involved in ripple formation. 

In conclusion, results from this feasibility study demonstrate that a 
wind tunnel of the carrousel design could be fabricated to operate in a space 
station environment and that experiments could be conducted which would yield 
significant results contributing to the understanding of aeolian processes. 
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FIGURE CAPTIONS 

1. Forces on a particle at the threshold of motion (here idealized as a 
sphere and resting on a ted like particles) include aerodynamic 
lift (L) > drag (D),and moment (M) , weight (W) and interparticle or 
cohesive force (Ip). The only force remaining to resist motion at 
zero gravity is the cohesive force Ip. Isolation of this force 
would permit much more accurate evaluation of its effect in aeolian 
proceiijses. 

2. Experimental data for threshold at the Earth’s surface show that the 
interparticle force Ip is the dominant resisting force for small 
particles and the weight is dominant for large particles, resulting 
in an "optimum" particle diameter for minimum threshold (from 
Iversen et al. , 1976). 

3. The difference between the two curves on the left side of the figure 
illustrates the approximate relative magnitude of the cohesive force 
compared to weight for smaller particles. 

4. Threshold friction speed versus particle diameter on Earth and the 
expectation of threshold under zero-g; large particles could be 
moved more easily, 

5. Swirl (dust-devil) threshold speed on Earth (one-gravity) compared 
to values expected at zero-g; as for the boundary-layer threshold, 
the swirl threshold speed would decrease for larger particles 

( right-hand-side) , 

6. Trajectories of windblown particles; the solid line is an 
experimental trajectory at the earth's surface, the broken lines are 
calculated (from White and Schultz, 1977). Zero-gravity experiments 
will have two major advantages; (1) the threshold speeds for larger 
particles are much lower, thus much lower wind and particle speeds 
are possible, and (2) elimination of the force due to particle 
weight means the aerodynamic forces can be measured (analysis of 
motion pictures) much more accurately than on Earth. Results would 
improve understanding of particle physics, having important appli- 
cations to problems involving blowing sand, dust, soil, and snow. 

7. Isometric representation of the carrousel drums. 

8. Photograph of prototype Carrousel Wind Tunnel (CWT). 

9. Dimensions of the inner and outer concentric drums of the prototype 

CWT. 

10. Close-up photograph of the hot-wire anemometer probe mechanism used 
to determine the flow properties of the prototype CWT. 

11. Wind velocity profile data at station 17.5 ram from tunnel side wall. 

12. Wind velocity profile data at station 39.5 ram from tunnel side wall. 


WinH velocity profile data at station 77.5 mm from tunnel side wall. 

Wind velocity profile data at station 126.5 mm from tunnel side wall. 

Comparison of theoretical wind velocity profiles at the four stations 
in CWT. The central-layer data are matched to the theoretical 
profile. The agreement of the four velocitjy profiles shows remarkably 
uniform wind flow across the test section. 

Wind speed ratios within the carrousel wind tunnel versus distance 
from the tunnel side wall. 

Particle trajectories in a zero-gravity vacuum corresponding to 
different Initial radial velocities. The 'coo ;dinate system is fixed 
to the particle launch point and rotates with the outer cylinder. In 
inertial space the trajectories are straight lines, but relative to 
an observer standing on the launch point, which is radially 
accelerating, the trajectories appear curved. 

Particle trajectory in zero-g, with one atmosphere of air, based on 
assumed drag characteristics. The outer cylinder is not rotating and 
the coordinate system is fixed to it at the launch point. The 
initial inward radial velocity of the 100 ym diameter particle is 
assumed to be 326.3 mm/sec., the friction speed of the outer 
cylinder. 

Induced wind velocity profile within the carrousel wind tunnel for a 
515 RPM difference in spin rates of the inner and outer drums. This 
profile was used to calculate the particle trajectories in Figs. ,18, 
20, and 21. 

The zero-gravity space carrousel particle trajectory (radial height 
vs. arc-length) is compared to the particle trajectory on the 
Earth’s surface for the same boundary layer velocity profile. 

Particle trajectory in the carrousel if it is on the Earth’s surface 
with the inner drum only rotating about a horizontal axis. The 
particle takes off at the lowest point on the outer drum surface. 
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